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Abstract 

Confinement of electrically charged test particles in the dilute plasma of mo- 
nopoles and pointlike magnetic dipoles is studied. We calculate the tension of the 
t^J- ■ string emerging between the infinitely separated test particles. The string tension is 

an increasing function of the dipole density provided other parameters of the plasma 
are fixed. The relevance of our results to confining gauge theories is discussed. 
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1 Introduction 
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We study the confining properties the plasma of the Abelian magnetic monopoles with a 
fraction of the magnetic dipoles in three dimensional Euclidean space-time. This kind of 
the plasma may appear in gauge theories such as the Georgi-Glashow model which pos- 
sesses the topoiogically stable classical solution called the 't Hooft-Polyakov monopole 
The dipoles are realized as the monopole-anti-monopole bound states. Since the long 
range gauge fields of the monopole are associated with an unbroken Abelian subgroup the 
long range properties of the monopoles and as a consequence, of the dipoles, are Abelian. 

The confining properties of the dilute Abelian monopole-antimonopole plasma are 
well known ||. The test particles with opposite electric charges are confined due to 
formation of a stringlike object between the charge and the anticharge. The string has a 
finite thickness of the order of the inverse Debye mass and a finite energy per the string 
length ("string tension"). Thus the potential between the test particles is linear at large 
distances 0. 

At a sufficiently high temperature the plasma of the Abelian monopoles undergoes 
the Berezinsky-Kosterlitz-Thouiess phase transition || f|, ||. The vacuum in the high 



temperature phase is filled with the neutral monopole-anti-monopole bound states ||, ||, [7j 
obeying non-zero magnetic dipole moments[|. The long range fields of the magnetic dipoles 
are weak and they can not induce a non-zero string tension between the electric charges 
separated by large distances. Thus the finite temperature phase transition separates the 
confining and deconfining phases of the model. 

Since the formation of the magnetic dipole states is supported by the increase of the 
temperature we expect to have a mixed plasma of the Abelian monopoles and dipoles at 
the confining side of the phase transition. In this paper we investigate the behavior of 
the electric charges in the mixed plasma. Below we show that the magnetic dipoles affect 
the string tension even despite of the inability of the pure magnetic dipole gas to confine 
electric charges. For simplicity the dipoles are considered to be pointlike. 

Note that in the finite temperature Georgi-Glashow model the charged bosons are 
supposed Jl2| to play an important role. In this paper we are interested in the pure 
monopole-dipole gas without any charged matter fields. The structure of the paper is as 
follows. The path integral formulation as well as simplest properties of the monopole and 
dipole gas are discussed in Section ^. The interaction of the electric charges are studied 
both numerically and analytically in Section ^. Our conclusions are summarized in the 
last Section. 



2 Gas of monopoles and magnetic dipoles 

The partition function of the monopole-dipole gas has the following form: 
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where ( m ((d) is the fugacity of the monopole (dipole) component of the gas, and 
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are the densities of the magnetic charges and magnetic dipole moments, respectively. We 
use the latin, a (greek, a) subscripts to denote the monopole (dipole) parameters. The 
magnetic charge of the a th individual monopole in units of the fundamental monopole 
charge, g m = 47r/g, is referred to as q a . In the dilute monopole gas the (anti)monopoles 
have unit magnetic (anti) charges, \q a \ = 1. 

The magnetic moment of the a th dipole is jl a = g m n a r a (no sum is implemented) , 
where r a is the "dipole size" and n a is the direction of the dipole magnetic moment. 

^^Note that the monopole binding is qualitatively similar to the formation of the instanton molecules 
in the high temperature phase of QCD suggested to be responsible for the chiral phase transition M. 
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For simplicity we consider a dilute gas of the point-like dipoles characterized only by 
a (fluctuating in general case) dipole moment and the space-time position. At the end 
of the paper we discuss the applicability of the results obtained in the pointlike dipole 
approximation to the real physical systems. 

In eq.(|l]) the integration over the dipole moments / d 3 p a is represented as the integra- 
tion over direction n a of the dipole moment and over the dipole size r a , weighted with a 
normalized distribution function F(r), / °° F(r) dr = 1. The function D(x) = 47TIXP 1 in 
eq.([l]) is the inverse Laplacian, and the subscript "reg" indicates that the monopole and 
dipole self-interaction terms are subtracted. 

The second line of eq.(|l]) can be rewritten as follows^]: 

J Dip ex v{-7^ J d 3 x (dip) 2 + i^wixa) + i^2(r a , d)cp(x a )} . (3) 

Substituting this equation into eq. ([!]), performing summation over the monopole charges 
q a and integrating over the dipole moment directions n a , and the monopole (dipole), x a ^ a ) 
positions, we get the following partition function for the gas of the monopoles and dipoles: 

Z = jDip exp{- J d 3 xC(x)} , (4) 

^ = 4^ - 2C ™ cos ^ - ^ J dr F < r > ^^wr 1 • < 5 > 

32tH J r \d,p\ 

where \a\ = \fd? . 

We study the magnetic monopole-dipole plasma in the weak coupling regime. The 
density of the monopoles and antimonopoles, p m , and the density of the dipoles, pa, can be 
calculated as the following expectation values taken in the statistical sum ([!]): p m =< M > 
and pd =< N >. In the weak coupling limit the leading order contributions to the 
(anti)monopole and dipole densities are proportional to the corresponding fugacities p], ||: 

Pm = 2(m , Pd = 47TCd • (6) 

The correlations of the fields ip in model (|J) are short ranged indicating that the 
photon gets a non-zero mass in the presence of the monopole plasma M. Indeed the 
expansion of lagrangian (|5|) in small fluctuations of the field ip gives: 

C = ^ U + y Cdl?) (dipf + ( m ip 2 + O(ip') , (7) 

where 

167T 



p? = g 2 m r 2 = — — / dr F(r) r 2 
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2 Here and below we omit inessential constant factors in front of the path integrals. 
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is the mean quadratic magnetic moment of a dipole. Thus the two-point correlation 
function < <f(0)if(x) > is exponentially suppressed at large distances as e~ mD ^, where 
m£> is the Debye mass in the monopole-dipole plasma: 

m% = e- 1 Ml, Ml = 3 -^. (9) 

9 

Here Mp is the Debye mass in the pure magnetic monopole plasma and e is the dielectric 
constant (permittivity) of the magnetic plasma, 

e = l + yCiiF- (10) 

The presence of the dipole states in the plasma leads to charge renormalisation, g 2 — > eg 2 
according to eqs.(|7]|J). In the next Section we show that the dipoles affect also the 
confining properties of the monopole plasma. 



3 Electric charges in magnetic monopole-dipole gas 

Let us consider two static infinitely heavy electrically charged test particles in the pure 
magnetic monopole plasma. The monopoles and anti-monopoles form a double layer at the 
sides of the minimal surface, spanned on the particles' trajectories M: at one side of the 
surface the magnetic density is positive while at another side the density is negative. The 
stringlike structure between the charges has a finite thickness of the order of the inverse 
Debye mass, and a non-zero energy density per unit area of the surface. Thus, 

the electric charges are confined in the magnetic monopole plasma. At large separations 
R between the test particles the confining potential is linear while at RMd <C 1 the 
potential becomes of the form const. R a , where a ~ 0.6, Ref. JlO . 



The behavior of the interparticle potential in the pure plasma of pointlike dipoles is 
different. The potential is linear at small interparticle separations, V(R) = const, pa f R, 
due to interaction of the dipole clouds gathered near the electric sources 0. At large 
distances R the potential is of the Coulomb type with the renormalized coupling g, g 2 —>■ 
eg 2 , where e is the dielectric constant given in eq.(]TT)|). 

Thus, the pure gas of pointlike dipoles gives rise to the nontrivial modification of 
the interparticle potential only at small separations R since the fields of an individual 
dipole is decreasing faster than the monopole field and we could expect only a short 
range modification of the interparticle potential. However, in the mixed magnetic dipole- 
monopole plasma the role of the dipoles at large distances becomes nontrivial. Physically, 
the monopole fraction of the plasma leads to the formation of the stringlike structure 
between electrically charged particles while the dipole fraction of the plasma interacts 
with the electric field inside the string. This is the way how the dipoles modify the linear 
term of the interparticle potential at large distances. Below we study the influence of 
the dipole fraction of the gas on the long-range potential between electrically charged 
particles. 
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Let us consider the contribution of the monopole-dipole gas into the quantum average 
of the Wilson loop operator, Wc, where contour C is the trajectory of the infinitely heavy 
test particle carrying the fundamental electric charge, g/2: 



rf[x 



W c = exp{i | J d 3 x (p{x) + (fr\x), B. 
The function rjc is defined as follows: 

rf{x) = it J d 2 a ij (y)e ijk d k D(x - y) 
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where the integration is taken over an arbitrary surface^] £c spanned on the contour C. 

Substituting eq.flTID into eq. QT|) and performing the transformations presented in Sec- 
tion ^ we derive the following representation for the contribution of the monopoles and 
dipoles to the Wilson loop: 
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where the expression in the square brackets has been normalized to zero at </? = and we 
have shifted the field ip — > ip — r] c . 

In the weak coupling regime the densities of the monopoles are small and the Wilson 
loop in the leading order is given by the classical contribution to eq. ( fTBD . The correspond- 
ing classical equation of motion is: 
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Below we consider the infinitely separated static charge and anticharge located at the 
points (y, z) = (±oo, 0). The x-coordinate is considered as a time coordinate. The string 
Sc in eq.(|T2|) is chosen to be flat and is defined by the equation z = 0. Thus, the classical 
solution of eq.flnp does not depend on x and y coordinates, (p c \ = ip c \(z). For the sake of 
simplicity we consider the case of the non-fluctuating dipole moments, G(r) = 5(r — r ) 
(below we use r instead of r ). 

The classical contribution to the Wilson loop ( |i~3D gives the area low, < Wc > gas 
const . exp{— Area cr c \}, where the classical string tension is 
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! Note that the Wilson loop does not depend on the shape of this surface. 
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and the ip c \ is the solution of the rescaled classical equation of motion (JT4|): 



<9fy?ci - sinv? c i - sdj: 



cos(fr|<9^ cl |) - 



2^5(0- (16) 



_|%?ci| \ h\d^ip ci \ t 

The dipole size in units of the Debye mass is denoted as 

h = rM D , (17) 

£ = z M D is the rescaled z-coordinate and 

o d 2ir( d 



•s 



is the ratio of the dipole density to the monopole density ("dipole fraction"), eq.@. 

In the absence of the dipoles, s — 0, the solution of eg. (|16D and the corresponding 
string tension are ||: 



^(0 = 4 sign® arctane-KI , a™ = -^H , (i 9) 
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respectively. In the presence of the dipoles the string tension is modified: 

a cl = a^H(s,h). (20) 

We study the function H below. 

In the limit of the small dipole size, h <C 1, eq.(|l6|) can be expanded in powers of h 
up to the second order and the classical string tension can be found analytically. The 
correction factor H is the square root of the dielectric constant flnp, 

H(s,h) = [l + ^(h 2 + 0(h 4 ))Y , (21) 

Obviously, this correction is due to the renormalisation of the coupling constant g 2 — > eg 2 . 

Unfortunately, eq.flTop can not be solved analytically for general values of the dipole 
fraction s and the dipole size h. Therefore we solve this equation and calculate the 
correction factor H numerically. The function H at various (fixed) dipole fractions s is 
shown in Figure |l|(a) by solid lines. The dashed lines indicate the behavior of correction 
factor H, eq.(pl[) valid at the small dipole sizes h. 

The behavior of function H can be understood as follows. The size of the flux of the 
electric field coming from the test particles is of the order of the inverse Debye mass. The 
dipoles with small size h do not feel the flux structure and the contribution of the dipole 
fraction of the gas to the string tension is solely due the polarization of the media of the 
magnetic dipoles by the electric field, eq. ([?!]). The correction factor if is a rising function 
of h since the increase of the dipole magnetic moment (size) leads to the enhancement of 
the mean polarization energy of the dipole in the external field. 
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Figure 1: (a) The correction factor H, eq.(^), to the string tension vs. the dipole size 
h for the different dipole fractions flTSP: s = 0.1, 0.2, 0.5; (b) the correction factor Q, 
eq. (p3[) , fs. the dipole fraction s for the dipole sizes h = 0.1, 1, 3. The numerical results 
are shown by the solid lines and the corresponding leading contributions at small h are 
shown by the dashed lines. 



At h ~ 1 the function H(h) starts to deviate from eq.(^Tj) since the dipoles begin to feel 
the width of the electric flux. At h « 2 the rise of the function H slows down, Figure |l|(a), 
and at large values of h the correction factor grows logarithmically, H ~ s log/i. To 
understand this effect we note that the pure dipole gas can not generate a finite correlation 
length || 11] contrary to the pure monopole gas. Therefore one can expect that the dipole 



fraction of the monopole-dipole gas has no essential influence on the Debye mass^. As a 
result, the derivative |<9,e<^ci| is of the order of unity at the center of the electric flux and 
the derivative vanishes out of the flux core as e~ c ^, C ~ 1. One can easily check that 
the dipole contribution to the string tension (the last term in eq. (|T5|) ) in the large h limit 
is equal to s at |£| <^ log/i and vanishes quickly out of this region. Since in this limit the 
other terms in eq.(^) do not grow with the enlargement of h, the string tension is given 
by the following formula, 



s 9 sJPr. 

a = 



(log/i + 0(l)), /i>l. (22) 



47T 

Another interesting question is the dependence of the string tension a on the dipole 
fraction s at fixed total number of the monopoles. The total (anti)monopole number in- 
cludes both the free ( ant i) monopoles and the constituent (anti)monopoles of the magnetic 

4 A numerical study of the profiles of the field (p c \(z) shows that this is indeed the case for all studied 
values of h and s. 
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dipole states. Thus, p tot = p m + 2p d = p m (1 + 2 s), where the dipole fraction s is defined 
in eq.([H|). If a part of the free monopoles transforms into the dipole states the string 
tension acquires positive a contribution due to the increase of the dipole states (the func- 
tion H grows) and a negative contribution due to the decrease of the monopole density 
(cr£i diminishes). The correction factor Q due to the monopole binding is defined by the 
following equation: 

a(s,h)=a tot Q(s,h), Q{s, h) = (1 + 2 S y 1/2 H(s, h) . (23) 

where a to t = 2g y/ptot / ^ is the string tension in the absence of the dipole states. Note 
that a to t remains constant during the binding process. 

In Figure |](b) we show the function Q vs. the dipole fraction s for a few fixed values 
of the dipole sizes h. In accordance with physical intuition the string tension (solid lines) 
decreases due to the monopole binding. Moreover, the larger the dipole size the smaller 
the decrease is. This is not an unexpected fact since the contribution of the dipole to the 
Wilson loop grows with the increase of the dipole magnetic moment (equivalently, "size" 
h of the dipole). 

Let us discuss the applicability of the pointlike dipole approximation to a physical 
system, for example, to the Georgi-Glashow model. The dipoles may be treated as point- 
like particles provided the typical distance between the plasma constituents is much larger 
than the dipole size r. The corresponding condition reads as follows: 

(p m + Pd) r 3 < 1 . (24) 

Using eqs. ( |T7|jrg| ) this equation can be rewritten as follows: 

(1 + s) h 3 n rn < 1 , n m = p m Mp 3 , (25) 

where n m is the number of the monopoles in the unit Debye volume. The string tension 
can be calculated in the classical approximation (|T5|JT6| ) only in the weak coupling limit 
which is equivalent to the condition n m ^> 1. One can easily check that the condition 
h 1 can be fulfilled in the Georgi-Glashow model and the correction factor (|2lD to the 
string tension can be arbitrarily large. The condition h ^> 1 can not be fulfilled and the 
result (|22D is not valid in the Georgi-Glashow model in the weak coupling regime. 

Conclusion and acknowledgments 

We have studied the confining properties of the gas of the Abelian magnetic monopoles 
and the pointlike dipoles. Both the increase of the dipole fraction at a fixed monopole 
density and the enlargement of the dipole magnetic moment amplify the coefficient in 
front of the linear potential ("string tension"). At the same time the monopole binding 
at a fixed total number of magnetic charges leads to the natural decrease of the string 
tension. 
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The described effects may have interesting physical applications since the monopole- 
dipole gases are realized in various gauge theories. For example, the monopole binding 
mechanism leads to the high temperature deconfining phase transition in the three di- 
mensional compact electrodynamics and in the Georgi-Glashow modef], Ref. j7|. An- 
other example is the electroweak model in which the formation of the Nambu monopole- 
antimonopole pairs has been observed [13 1 in the low temperature (Higgs) phase. In this 



model the string tension for the spatial Wilson loops is non-zero in the high temperature 
(symmetric) phase in which the monopoles form a gas. At low temperatures the string 
tension vanishes |T3[ . 



Author was partially supported by RFBR 99-01230a and Monbushu grants, and CRDF 
award RP1-2103. 
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